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ABSTRACT 

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A -with 
the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, 
and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant 
signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play 
an essential role for the separation between 3C 66 A and the nearby radio galaxy 3C 66B, which is 
at a distance of only 6'. The derived integral flux above 100 GeV is 8.3% of Crab Nebula flux and 
the energy spectrum is reproduced by a power law of photon index 3.64 ± 0.39stat i 0.25sys. Within 
errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected 
for absorption by the extragalactic background light, we only find small differences between the four 
models that we applied, and constrain the redshift of the blazar to z < 0.68. 

Subject headings: BL Lacertae objects: individual (3C 66A) - galaxies: active ~ gamma rays: galaxies 
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1. INTRODUCTION 

Blazars make up the majority of extragalactic sources 
of very high energy (VHE; E > 100 GeV) gamma rays. 
They are a subset of active galactic nuclei (AGNs), 
and consist of BL Lacertae (BL Lac) objects and flat- 
spectrum radio-loud quasars. The general framework to 
explain the gamma ray emission is that they are pro- 
duced by charged particles which are accelerated in a 
relativistic jet. These jets are powered by gas accretion 
into a central supermassive black hole and are perpen- 
dicular to the accretion disc. When the jet is directed to 
us, the energy and flux of gamm a-rays are boosted by the 
relativ istic beaming effect ('e.g.. lBlandford fc Reeslll978l : 
lUrrv fc Padovani 1995). 

Generally, the spectral energy distribution of AGNs 
can be described by two broad bumps. The lower ener- 
getic bump, at frequencies from radio to X-rays, is at- 
tributed to synchrotron emission from nonthermal rel- 
ativistic electrons in the jet. The other bump, cov- 
ering the X-ray to gamma-ray bands, could either be 
due to inverse Compton scattering of seed photons 
by the electrons (leptonic rnodel, e.g., [M araschi et al. 
1992t iDermer fc Schlickeiseii 119931: iBloom fc Marschen 



1996 ; lKrawczvnskill2004 ) or d ue to hadronic interactions 
(see, e.g., iMannhejml [19931; iMiicke fc Protherod 120011: 
iMiicke et al . 2003). 

3C 66A was class ified as a BL Lac object by 
iMaccagni et al.l (|1987l ). based on its significant optical 
and X-ray variability. The synchrotron p eak of this 
sourc e is located between 10^^ and 10^® Hz ijPerri et al.l 
[200l . therefore 3C 66A can also be classified as an 
intermediate- frequency peaked BL Lac object (IBL). The 
redshift of 3C 66A w as determined to be 2 = 0.444 by 
indep endent authors ([Miller et al.lll978l iLanzetta et al.l 
[19931) . However, their measurements are based on the de- 
tection of one single line. Another observatio n of 3C 66A 
at a d ifferent spectral range was reported bv [Finke et al.l 
([20081 ). but no spectral feature was found, and a lower 
limit of the redshift was derived to be 0.096. For the 
marginally resolved host galaxy ([Wurtz et al.l[1996i ). a 
redshift of 0.321 was found. Recently, through the inves- 
tigation of the Large Area Telescope (LAT), on board the 
Fermi Gamma-ray Space Telescope (Fermi) satellite and 
VHE gamma-ray observations, upper li mits for the red 
shift of 3C 66A were derived; z = 0.44 ([Prandini et al 
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201c , 2 a confidence level) and z = 0.58 ([Yang fc Wang 



Several gamma-ray observations of 3C 66A were per- 
formed since the 1990s. With the EGRET satellite, a 
GeV gamma-ray e mission (3EG J0 222+4253) was asso- 
ciated with 3C 66A ([Hartman et al.l fl999). However, due 
to the large EGRET pomt-spre&d function (PSF), an in- 
fluence by th e nearby pulsar PSR J0218-I-4232 could not 
be excluded ([Kuiper et al.|[2000l ). The Crimean Astro- 
physical Observatory claimed detections of 3C 66A above 
900 GeV with an integral flux of (3 ± 1) x 10"" cm^^g-i 
(fStcp an van et al.l (2002,) . Later observations by HEGRA 
and Whipple reported upper limits of f (> 6 30 GeV) < 
1.42 X 10-"cm-2s-i ([Aharonian et al.''2000') an d F{> 
350 GeV) < 0.59 x 10"" cm'^s"! (Horan ct"aL| [200l , 
respectively. Additionally, the STACEE observation 
found a hint of signals at a 2.2 significance level and 
derived upper limits of < 1.0 x lO^^^cm^^s^^ and 



< 1.8 X 10~"cm~^s~ ^ for threshol ds of 147 GeV and 
200 GeV, respectively ([Bramel et al.i r2005). 

Recent VERITAS observations of 3C 66A taken from 
2007 September to 2008 January and from 2008 Septem- 
ber to 2008 November, for a total of 32.8 hr, re sulted 
in a detection in VHE gamma rays ([Acciari etan[2009). 
The energy spectrum was derived with a photon index 
of F = 4.1 ± 0.4stat ± 0.6sys. The integral flux of the 
VERITAS observations above 200 GeV is (1.3 ± 0.1) x 
10-" cm-2 s-i (6 % of the Crab Nebula flux). 

3C 66A has been monitored by Fermi/LAT since 2008 
August, covering the latt er part of the VER ITAS ob- 
servation. According to lAbdo et"an ([20091) . who re- 
ported the flrst 5.5 months of Fermi/LAT observations 
of 3C 66A, the blazar showed a significant fiux variability 
(a factor of 5-6 between the highest and lowest fluxes). 
The derived energy spectrum with the photon index of 
F = 1.98 above 1 GeV, in combination with the VERI- 
TAS spectrum, indicates that the spectrum must soften 
above 100 GeV. 

MAGIC observed the sky region around 3C 66A from 
2007 August to December, obtainin g a total exposure 
time after data quality cuts of 45.3 hr ([Aliu et a l. 2009b). 
These data revealed a significant VHE gamma-ray sig- 
nal centered at 2'^23'"12", 43°0'7". This excess (named 
MAGIC J0223-I-430) coincides within uncertainties with 
the position of a nearby, Fa naroff-Riley-I (FR I) type 
galaxy 3C 66B (z = 0.0215: [StuIFeTal] [19751) . StiU, 
judging from the skyplot alone, the probability of the 
emission to originate from 3C 66A is 14.6 %. The en- 
ergy spectrum of MAGIC J0223-I-430 was reproduced 
by a single power law with the index of F = 3.1 ± 
0.3. The integral flux above 150 GeV corresponded to 
(7.3 ± 1.5) X 10"^^ c m-^ s"^ (2.2 % of the C rab Nebula 
flux). According to [Tavecchio fc Ghisellin il p008), the 
radio galaxy is also a plausible source of VHE gamma- 
ray radiation. Also, the r ecent MAGIC detection of IC 
310 ([Mariotti et al. |[2010( ). a radio galaxy at a very sim- 
ilar redshift (z = 0.0189) indicates that 3C 66B might 
be feasible to explain all or part of the MAGIC detection 
from 2007. 

2. OBSERVATIONS 

From mid 2009 August, 3C 66A went into an optical 
high state which was reported by the Tuorla blazar mon- 
itoring prograrrF^. This outburst triggered new MAGIC 
observations. The optical flux in the R band reached 
a maximum level of ~12mJy in 2010 January, while 
the baseline flux in the historical data of the source is 
'-~^6 mjy. 

The observations were carried out with the MAGIC 
telescopes located on the Canary Island of La Palma 
(28?8 N, 17?8 W, 2220ma.s.l.). The two 17m diame- 
ter telescopes use the atmospheric Cherenkov imaging 
technique and allow for measurements at a threshold as 
low as 50 GeV in normal trigger mode. 

We observed the blazar 3C 66A in several time slots 
between 2009 September and 2010 January. However, 
the sky imaging CCD cameras that are used to cross- 
check the telescope pointing ("starguider cameras") only 
became fully applicable to stereo observations in early 
December. To allow for a high-confldence directional 

http:/ /users. utu.fi/kani/lm/indox. html 



Observations of the Blazar 3C 66A with the MAGIC Telescopes in Stereoscopic Mode 



3 



statement on the arcminute scale, we therefore only used 
data taken after these upgrades, which were 5.6 hr in 
total. Furthermore, we had to discard data with low 
event rates, affected by the exceptionally bad weather 
conditions in that winter. Finally, we had 2.3 hr of good 
quality data left after all quality cuts. They were taken 
on six days between 2009 December 5 and 2010 January 
18, partly under low-intensity moon light conditions. 

The data were taken using the f alse source tracking 
(wobble) method (|Fomin et al.lll99^ , in which the point- 
ing direction alternates every 20 minutes between two 
positions, offset by ±0?4 in RA from the source. These 
wobble positions were chosen with respect to 3C 66A, but 
the small distance to 3C 66B (0?01) allows equal judg- 
ment for both sources. The data were taken at zenith 
angles between 13° and 35°. 

3. DATA ANALYSIS 

For the analysis, only stereoscopic events triggered by 
both MAGIC telescopes were use d. They were analyze d 
in the MARS analysis framework (|Moraieio et al.ll2009l ). 
taking advantag e both of the adv anced single-telescope 
algorithms (e.g.. lAliu et al.ir2009al ) and newly developed 
stereoscopic analysis routines. These routines are at 
present still subject to some minor improvements and 
will be discussed in more detail in a separate paper still 
in preparation, but are shortly outlined in the following. 

Combining monoscopic and stereoscopic strategies, the 
direction of gamma rays is calculated for each tele- 
scope separately, u sing the random forest technique 
(| Albert et al.ll2008bl ). and later combined with the pro- 
jected crossing point of the image axes, with a weight 
depending on the angle between the two shower images. 
Requiring a certain level of agreement between the differ- 
ent estimates furthermore improves the resolution, and 
also helps to reject the (less focused) hadron showers. 
Similarly, an energy estimator is determined from look- 
up tables for each telescope separately, and later com- 
bined to a common estimated energy. 

The skymap generation, which is particularly impor- 
tant for the analysis of data from the 3C 66A/B region, 
follows a two-step algorithm. The first step is to gen- 
erate an exposure model for the field of view in camera 
coordinates, for the quality cuts that were applied in the 
analysis. This is done by joining the distributions of 
photon-like events from the two wobble positions, taking 
advantage of the fact that the source, in relative camera 
coordinates, is on opposite sides for both wobble sets. 

The second step is the calculation of an expected back- 
ground event distribution in celestial coordinates, and its 
comparison to the actual event distribution. Before that 
comparison, a smearing with a Gaussian kernel is ap- 
plied. The significance s are calculated following Equa- 
tion (17) of Li fc Mai (|1983l ). taking into account the 
higher precision of the background estimation implied 
by the above modeling. 

The performance of the analysis software was opti- 
mized and checked with contemporaneous Crab Nebula 
data and MC. The Crab Nebula spectrum could be an- 
alyzed down to about 50GeV, fully covering the range 
of the spectrum presented in the next paragraph. The 
achieved angular resolution, defined as the ct of a two- 
dimensional Gaussian function, is around 0?1 at 100 GeV 
and approaching 0?065 at higher energies. This a defines 



I 3C 66A/B Significance Map [ 




2.42 2.4 2.38 2.36 2.34 

RA [h] 

Fig. 1. — MAGIC significance skymap of the region around 
3C 66A/B for events with energies above 100 GeV. 

the radius in which 39 % of all photons of a point source 
are contained. The systematic uncertainty on the direc- 
tion reconstruction is a product of the telescope pointing 
uncertainty and possible biases that occur in the recon- 
struction algorithms. The latter can be caused by irreg- 
ularities in the shower images, such as missing camera 
pixels, inhomogeneous noise from stars in the field of 
view, or imperfections in the data acquisition electron- 
ics. Both the total pointing deviation and the telescope 
pointing precis ion of MAGIC were always m onitored over 
the years (Br etz et afl 120091 : lAleksic et a l. 2010), and 
along with studies of contemporary stereo data of known 
direction lead to an estimate of the maximal systematic 
stereoscopic pointing uncertainty of 0?025. 

We also used the publicly accessible f ermi/LAT datsF^ 
to investigate the status of the source in the GeV energy 
range during the MAGIC observation period. The Fermi 
data were analyzed using the public software package 
LAT Science Tools v9.15.2, including the Instrument Re- 
sponse File P6_V3_DIFFUSE, and galactic, extragalactic 
and instrumental background models. 

4. RESULTS 

Figure [I] shows a skymap of the observed region above 
100 GeV. The significance of the excess at the location 
of 3C 66A is 6.4(7. We cross-checked the detection also 
by investigating the distribution of squared angular dis- 
tances {9'^) between photon directions and the assumed 
source position. The expected background is extracted 
from corresponding 9'^ plots done with respect to other 
sky positions at similar distance from the pointing di- 
rection. Comparing the data with this expectation we 
find a significance of 5.2 cr (see Figure [2]). The difference 
in significance can be attributed to the different integra- 
tion procedure of signal and background in the skymap, 
which generally leads to a slightly better background es- 
timation and therefore a higher significance. 

We also analyzed the data taken with and without 
moon light separately to find possible effects from the 
higher thresholds of individual camera pixels. However, 
we could not find a clear tendency beyond the statisti- 
cal errors and thus decided to use all the data for the 
analysis. 

http;//fcrmi. gsfc.nasa.gov/ 
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Fig. 2. — Distribution of squared angular distances between pho- 
ton directions and the position of 3C 66A (6^) for events with 
energies above 100 GeV. The OFF data are taken from three po- 
sitions that are symmetrical with respect to the telescope pointing 
directions. 

Unlike in the 2007 observations of this sky region, the 
emission peak this time is clearly on top of 3C 66A. The 
fitted center of gravity of the excess (small black square 
in Figure H]) is at a distance of 0?010 ± 0?023 (stat.) ± 
0?025 (sys.) from 3C 66A, and 0?108 ± 0? 023 (stat.) ± 
0?025 (sys.) from 3C 66B. While being compatible with 
the former, the statistical rejection power for the emis- 
sion to emerge from the radio galaxy 3C 66B corresponds 
to 4.6 standard deviations. Even considering the unlikely 
case of a systematic offset exactly toward the blazar, the 
rejection significance of 3C 66B is at least 3.6 cr. These 
numbers were confirmed by a second analysis with inde- 
pendent data quality selection and cut optimization pro- 
cedures. The same result is found even when the photon 
direction is taken only from the projected crossing point 
of the two shower axes. We therefore conclude that the 
signal we see this time emerges from the blazar 3C 66A. 

It shall be mentioned that this result is a clear merit 
of the angular resolution and background rejection of 
the new stereoscopic system. In fact, if we compare the 
above stereo directional reconstruction algorithm to the 
MAGIC-I algorithm alone, we find basically the same 
result, but the statistical error of the fitted source posi- 
tion increases roughly by a factor of two. Consequently, 
the rejection significance of 3C 66B would be less than 2 
standard deviations, and the total detection significance 
would be below 5 standard deviations. 

The energy spectrum of 3C 66A was derived using four 
different unfolding algorithms (Albert et al. 2007) which 
correct for efficiency, smearing and biasing effects in the 
energy response of the detector. The most conservative 
of these methods is the so-called forward unfolding, in 
which essentially a spectral shape is assumed a priori, 
and its parameters are adjusted by iteratively folding 
the assumed spectrum with the response function until 
the predicted distribution of estimated energies matches 
optimally the actually measured distribution. With all 
unfolding methods, we found that the data are well com- 
patible with a power law of the form 



dF 

dE ^ 



E 



200 GeV 



-r 



(1) 



with a photon index F 
flux constant at 200 GeV of ^^200 



3.64±0.39stat±0.25sys and a 
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Fig. 3. — Observed and EBL-corrected (de-absorbed) differential 
energy spectra of 3C 66A in the period of 2009 December and 2010 
January. The light-shaded area indicates the 1 a range of the ob- 
served power law spectrum gained by forward unfolding (see the 
text) , the crosses are from the unfolding after Tikhonov &: Arseni^ 
II1979I) for comparison. The dark-shaded area is the spread of the 
de-absorbed, mean flux values obtained by the four applied EBL 
models, assuming the redshift of z=0.444 . The VERITAS (ob- 
served) spectrum after IXcciari et al.l 1I2OO9I') is shown for compari- 
son. 

corresponds to (4.5 ± 1.1) x 10"" cm^^ g-i (8.3 % Crab 
Nebula flux) . Here, the parameters and statistical errors 
are taken from the forward unfolding, while the system- 
atic errors reflect the variations among the other unfold- 
ing algorithm s, plus several s tandard uncertainties dis- 
cussed in Albe rt et al.l ()2008a[ ). The systematic flux un- 
certainties add up to 36 % in total. Figure [3] displays the 
function we fitted through forward unfolding, and spec- 
tral points derived using th e Tikhonov unfolding method 
(jTikhonov fc Arsenin|[l979l ). 

Due to the shortness of our observation, we cannot 
discuss flux variability with these data. However, com- 
paring the flux to the one from our previous observation 
of the 3C 66A/B region confirms the VERITAS report 
of 3C 66A being a variable source in general. 

We also analyzed the Fermi data from the same time 
period. The fiux variability we found in a week-to-week 
light curve is not significant. Given the statistical un- 
certainties of the light curve, we would be sensitive on 
3(7 level to fiux variations of 60% or greater, and con- 
clude the variability in the days we observed must be less 
than that. The averaged flux above 200 MeV is roughly 
comparable to the averaged fiux over the first 5.5 months 
()Abdo et al.ll2009f ). and lower than that seen in 2008 Oc- 
tober, when a strong TeV flare was observed by VERI- 
TAS. A single power- law model can reproduce the source 
spectrum, a nd the photon ind ex is compatible with the 
one found in lAbdo et al.l (|2009[ ). indicating no signiflcant 
change in the overall spectral shape. 

5. DISCUSSIONS AND CONCLUSIONS 

MAGIC observed the 3C 66A/B region in 2009 De- 
cember and 2010 January, during an optical active state 
of 3C 66A and detected a clear VHE gamma-ray sig- 
nal. The excess coincides with the position of 3C 66A, 
and we rule out the emission to come from 3C 66B at a 
confidence level of 3.6 cr. This detection does not contra- 
dict the earlier MAGIC detection, though, which favored 
3C 66B as the VHE source. On the one hand, because 
the observation time of 2.3 hr would be too short to de- 
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TABLE 1 
EBL Corrected Indices 



Model Tint 

Franceschini et al. (20081 2.57 ± 0.68 

Gilmorc et al. (20097^ 2.61 ± 0.67 

'DommKucz ct al. (2010) 2.59 ± 0.68 

Kneiske fc Dole (20loT 2.37 ±0.70 



tect the VHE emission of 3C 66B, if on a similar flux level 
as in 2007, and on the other hand, because its flux may 
be even lower than before. In fact, 3C 66A might have to 
be in a low flux state in order not to outshine the compa- 
rably weak emission from 3C 66B at this close distance 
of about 1 cr of the PSF of the MAGIC telescopes. 

The obtained energy spectrum is softer than in the 
previous MAGIC detection (T = 3.10 ± 0.31stat ± O-'^sys) 
and compatible with the VERITAS spectrum of 3C 66A. 
Compared to VERITAS, the MAGIC measurement has 
a lower threshold and the spectrum is extending to well 
below 100 GeV. The flux level of 8.3 % Crab Nebula flux 
is similar to the one reported by VERITAS (6%), and 
significantly higher than in the previous MAGIC obser- 
vation (2.2%). 

The VHE photons produced at the source can be 
absorbed in the intergalactic space by pair production 
with the low energy (UV to infra red) photons of ex- 
tragalactic backgrouii d light (EBL; iStecker et al.l 119921 : 
iHauser &: Dwekl 720011 ). The amount of absorption de- 
pends on the energy and redshift, and can be cor- 
rected for in the data, assuming a certain modeling 
of the EBL density. Such a de-absorbed spectrum 
can be regarded as the spectrum we would measure 
if there were no EBL. To derive a de- absorbed spec- 
trum, we tested several state of the art EBL models, 
na melv. iFrancesc hini ct al . (20 d 8l). the fi d ucial model 
in iGilmore et al.l (120091) ■ iKneiske fc Dold (fMot ). and 
iDommguez et all (|20100 . The EBL corrections were ap- 
plied in the spectrum unfolding procedure (see above), 
using the full covariance matrix to correctly calculate the 
errors. The spread of the differential, de-absorbed flux 
spectra, obtained with the four models and assuming the 
redshift of z = 0.444, is shown as the dark shaded area 
in Figure [21 The de-absorbed photon indices for the four 
EBL modelings are listed in Table [TJ The differences be- 
tween the de-absorbed spectra are very smal l, although 
the one corrected after IKneiske fc Dold (|2010l ) is slightly 
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harder than the others. This also reflects the fact that 
also the predicted EBL shapes and densities are very sim- 
ilar in the first t hree m odels, but the overall density in 
IKneiske fc Dold (|2010( ) is somewhat higher. 

From most VHE emission models, the de-absorbed 
spectrum is expected not to be concave, i.e., rising to- 
ward higher energies. This can be tested both by com- 
paring the points of our own spectrum, but also by a 
comparison with the Fermi photon index (1.98). The 
fact that we find our spectrum neither significantly con- 
cave nor harder than in Fermi suggests that the as- 
sumed redshift of z = 0.444 does not contradict our 
observations. In fact, we investigated the plausibil- 
ity of the redshift, assuming that the intrinsic spec- 
trum is not expected to be exponentially rising, and 
thus have a pileup, at highest energies. This common 
method was previously used and described, for exa mple, 
in Maz in fc Goebel (2007.] : .Mazin fc Raud (pOOl . Us- 
ing the IFranceschini et al.l ()2008D model and the likeli- 
hood ratio test between the " power law" and " power law 
-I- pile-up" hypotheses, as described in the reference, we 
derive an upper limit on the redshift of z < 0.68. 

The results derived in this paper demonstrate the ad- 
vantages of the MAGIC stereoscopic system. Further 
MAGIC and other gamma-ray observations of this re- 
gion can provide interesting information about the IBL 
type BL Lac object 3C 66A, and, during low fiux periods 
of that, also the FRI type galaxy 3C 66B. 
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